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^ ' Abstract 

<D ■ 

. In this paper we consider a class of planar autonomous systems having an isolated limit cycle Xq of smallest 

. period T > such that the associated linearized system around it has only one characteristic multiplier 

with absolute value 1. We consider two functions, defined by means of the eigenfunctions of the adjoint of 
the linearized system, and we formulate conditions in terms of them in order to have the existence of two 
geometrically distinct families of T— periodic solutions of the autonomous system when it is perturbed by 
nonsmooth T— periodic nonlinear terms of small amplitude. We also show the convergence of these periodic 
solutions to xq as the perturbation disappears and we provide an estimation of the rate of convergence. The 
employed methods are mainly based on the theory of topological degree and its properties that allow less 
regularity on the data than that required by the approach, commonly employed in the existing literature on 
this subject, based on various versions of the implicit function theorem. 

00 
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' 1. Introduction 

Loud in [23] provided conditions under which the perturbed system of ordinary differential equations 

•l-H 

j_j ■ x = ip(x)+e(f)(t,x,e), (1) 



(N 
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u 



c3 



where 

t/> G C 2 (M™,R n ), (j) G C\R x R" x [0,1], M") (2) 

and 4> is T-periodic with respect to time, has, for sufficiently small e > 0, a T-periodic solution which tends to 
a T-periodic limit cycle xq of the unperturbed system 

x = ip(x) (3) 

as £ — > 0. The limit cycle xq satisfies the property that the linearized system 

y = i>'(x (t))y (4) 
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has only one characteristic multiplier with absolute value 1. Here and in the following by C*(R m , R n ) we denote 
the vector space of all continuous functions acting from R m to R™ having i-th. continuous derivatives. The main 
tool employed by Loud is the following, so-called bifurcation, function 

/o(0) = / T (zo(r),0(T-0,xo(T),O))dr, (5) 
Jo 

where z is a T-periodic solution of the adjoint system of (4) 

z = -W(x (t)))'z, (6) 

here A* denotes the transpose of the matrix A. Specifically, Lemma 2 in [23] states that in order that system 
(1) has a T-periodic solution x e such that 

x e (t - O ) -> xo(t) as e -> (7) 

it is necessary that 0o G K be a zero of the equation 

MO) = 0. (8) 

If (8) is satisfied for some = O and /o(0o) ^ 0, i- e - $o is simple, then by ([23], Theorem 1) for all sufficiently 
small e > system (1) possesses a T-periodic solution x e satisfying 

\\x £ (t-0 o )-x o (t)\\<sM, (9) 

where M > is a constant. These results are also consequences of general results stated by Malkin in [26]. 

The function / has been widely employed to treat different problems concerning periodic solutions of system 
(1) with e > small. We quote in the sequel some papers from the relevant bibliography devoted to this subject. 
In [23] Loud also considered the case when (8) is identically satisfied, i.e. /o(0) = for any e [0, T], to treat this 
case he introduced a new function which plays the role of fo and he showed that if this function has a simple zero 
0o then there exists a family of T-pcriodic solutions to (1) satisfying (7) (see also [22]). Moreover in [23] it is also 
considered the case when O is not a simple zero of /o, and the problem of the existence of T-periodic solutions 
to (1) is associated with the problem of the existence of roots of a certain quadratic equation. The case when 
the limit cycle xo of system (3) is not isolated, in particular, when the unperturbed system is Hamiltonian, and 
the case when the characteristic multiplier of system (4) is not simple have been considered by many authors. 
If system (3) is not necessary autonomous and it has a multi-parameterized family of T-periodic solutions, then 
existence of T-periodic solutions of the perturbed system satisfying (15) was proved by Malkin [26]. Mclnikov 
[29] treated the case when the limit cycle is not isolated and the limit cycles near xq are of different periods 
(see also Loud [24] and Kac [14]) and he showed that the simple zeros of suitably defined bifurcation functions 
/to, rn m,n e N, called Melnikov subharmonic functions, generate periodic solutions in a neighborhood of Xo 
whose periods arc in m : n ratio with respect to the periods of the perturbation term. Finally, Rhouma and 
Chiconc [34] have considered the case when 1 is not a simple multiplier of the linearized system, to deal with 
the problem of existence of T-periodic solutions they introduced a new two variables bifurcation function / 
whose simple zeros determine families of T-periodic solutions satisfying (7) . 

These theoretical results have been then developed in different directions: Hausrath and Manasevich [10], 
(see also [11]), found a class of T-periodic perturbations <j> for which the subharmonic Mclnikov function /^i 
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has at least two simple zeros, obtaining the existence of at least two families of T-periodic solutions to (1) 
satisfying (7). Makarenkov in [27] provided useful formulas to calculate simple zeros of Malkin's bifurcation 
function in case when the function <f> is sinusoidal in time. Tkhai [37] and Lazer [20] developed Malkin's 
and Mclnikov's approaches respectively to study the existence of periodic solutions to (1) satisfying (7) and 
possessing some additional symmetry properties that represent relevant features in the applications. Farkas in 
[13] investigated the existence of the so-called D-periodic solutions to (1) which arc not necessarily periodic 
but having periodic derivative. Greenspan and Holmes in [8] and Guckenheimer and Holmes in [7] applied the 
method of subharmonic Melnikov's functions to a variety of practical problems, a number of applications of 
Malkin's bifurcation function can be found in the book of Blckhman [1]. 

In all the previous papers, to show the existence of T-periodic solutions for e > small, several formulations 
of the implicit function theorem have been employed. Therefore, condition (2) is the common assumption 
of these papers (sometimes it is even required more regularity on ip and (f>). The persistence of the limit 
cycle x under less restrictive regularity assumptions than (2) is studied only for the cases when system (3) is 
linear, in this case the modified averaging methods developed by Mitropol'sksii [30] and Samoylenko [35] can be 
applied as well as the coincidence degree theory introduced by Mawhin, see, for instance, ([28], Theorem IV. 13); 
Hamiltonian, see M. Henrard and F. Zanolin [12]; or piecewise differentiable, see Kolovskii [17] and Steinberg 



In the present paper we assume that the linearized system (4) has only one characteristic multiplier equal 
to 1 and 



where z\ is an eigenfunction of system (6) corresponding to the characteristic multiplier /?» ^ 1, we give 
conditions in Theorem 3 for the existence of T-periodic solutions to (1) satisfying (7). Although, as we have 
mentioned before, in many papers it was proved the existence of two or more families of T-periodic solutions to 
(1) converging to xo in the sense of (7), it was not guaranteed that these families do not coincide geometrically, 
namely if one is just a shift in time of the other. In this paper our results ensure the existence of at least 
two geometrically distinct families of T-periodic solutions to (1) satisfying (7). Moreover, since property (9) 
is a consequence of the application of the implicit function theorem it is not anymore guaranteed under our 
conditions (10). However, we will show in Theorem 1 that under conditions (10) the following property holds 

eMi|/i(0 o ,t)| < \\x £ {t-e e {t)) ~Xo{t)\\ +o{e) < eM 2 \fi(6 ,t)\ for any e G (0, e ) and any t E [0,T], (11) 

where < Mi < M 2 and e (t) — > 9q as e — ► uniformly with respect to t E [0,T]. The introduction of the 
function /i, as shown by (11) and Corollaries 1 and 2 of this paper, gives a new qualitative information about 
the convergence (7) with respect to (9) and it is a contribution to the problem posed by Hale and Tboas in [9] 
concerning the behavior of the periodic solutions of a second order periodically perturbed autonomous system 
when the perturbation disappears. We would like also to remark, that Loud in [23] provided a precise information 
about the way of convergence of x e to x by means of the representation x £ {t) = xo(t + 9 ) + ey(t + 9q) + o(e), 



[36]. 



tp E C^R^K 2 ), (j) E C(R x R 2 x [0, 1],M 2 ). 



(10) 



By combining the function / with the analogously defined function 




3 



where the function y is a suitably chosen solution of system (60) of this paper with £ = xo(0), sec ([23], formulas 
1.3 and 2.11). 

In order to prove the existence of T-periodic solutions to (1) satisfying (7) under assumptions (10) we make 
use of the topological degree theory. Specifically, for e > 0, we consider the integral operator G e : C([0, T], R 2 ) — > 
C([0,T],R 2 ) given by (G E x)(t) = x{T) + f* i/)(x(r))dT + £ f* <p(r, x(t) , s)dr, t G [0,T], here C([0,T],R 2 ) is the 
Banach space of all the continuous functions defined on [0, T] with values in R 2 equipped with the sup-norm. We 
also consider the Leray-Schauder degree d(I — G e , Wu), see Brown ([4], §9), of the compact vector field I — G e 
with respect to the open set Wu = {x e C([0, T], R 2 ) : x(t) G U, for any t G [0, T]} , where U is an open set of 
R 2 . We will provide conditions in terms of the functions /o and /i ensuring that d(I— G s , Wu ) ^ d{I — G e ,Wu e ) 
for all e > sufficiently small, where Uq is the interior of the limit cycle xo and U £ C Uq (or Uq C U £ ) is a 
suitably defined family of sets such that U e —* Uq as e — > 0. To do this we use a result by Capietto, Mawhin and 
Zanolin ([5], Corollary 1) which, under our assumptions, states that d(I — G , Wu E ) — 1 for e > sufficiently 
small. Then in Theorem 2, by means of a result due to Kamcnskii, Makarenkov and Nistri ([15], Theorem 2), 
we conclude that there exists a continuous vector field F : R 2 — > R 2 with 

F(x o (0)) = M6)±o(6) + fi(0,0)y 1 (0), for any G [0, T], (12) 

here y\ denotes the eigenfunction of (4), corresponding to the characteristic multiplier p ^ 1, such that for all 
e > sufficiently small we have that d(I — G e , Wu a ) — ds{F, Uq), where ds(F, Uq) is the Brouwer degree of F 
on Uo, see e.g. Brown ([4], §8). In our case the integer ds{F, Uq) can be easily calculated since it is equal to the 
Poincare index of x with respect to the vector field F multiplied by +1 or —1 according with the orientation 
of xo, see Lefschetz ([21], Ch. IX, §4). Furthermore, as it was observed by Bobylev and Krasnoselskii in [2], 
for any small neighborhood Bs(dWu ) of the boundary dWu of Wu 01 we have that d(I — Go, Bs(dWu )) = 
and so one cannot directly apply Leray-Schauder fixed point theorem for studying the existence of T-periodic 
solutions to (1) satisfying (7). 

The paper is organized as follows. Theorem 1 of Section 2 states property (11) for the T-periodic solutions of 
system (1), in Theorem 2 we prove the coincidence degree formula d(I — G e , Wjj a ) = ds{F, Uq) for e > small. 
Finally, we give the main result of the paper: Theorem 3 which states the existence of at least two geometrically 
distinct families of T-periodic solutions to (1) satisfying (7). In Section 3 we provide an example which shows 
how formula (12) can be used for the practical calculation of ds{F, Uq). In fact, under quite general conditions 
on f and /i, we show that if </>(t, = -<p(t + T/2, £) for any t G [0, T] and £ e M 2 , then d B (F, Uq) G {0, 2}. 
Finally, we outline some methods for calculating the eigenfunctions j/i, zq and z\. 

2. Main results. 

Through the paper we assume the following condition: 

(A )— system (3) has a limit cycle x with smallest period T > and the linearized system (4) has only one 
characteristic multiplier equal to 1. 

In what follows we provide the notations that we will use in the proofs of the results of this Section. 

By yi we denote the eigenfunction of (4) corresponding to the characteristic multiplier p ^ 1 (clearly xq is 
the eigenfunction of (4) corresponding to the characteristic multiplier 1). Moreover, zq and z\ will denote the 
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eigenfunctions of (6) corresponding to the characteristic multipliers 1 and p* ^ 1 respectively. (01,02) is the 
matrix whose columns are the vectors 01,02 G R 2 , (ai,o 2 )* denotes the transpose of (01,02), fl(-,to,£) is the 
solution of system (3) satisfying x(t ) = £, 0^(-,fo,0 is the derivative of Q(-,t 0} £,) with respect to the third 
variable, U is the interior of the limit cycle x of system (3), dU is the boundary of U , [v]i, i = 1, 2, is the i-th 
component of vector v G R 2 . a || will indicate that the vectors a, b G R 2 are parallel, a 1 - denotes the vector 
a G R 2 rotated of ir/2 clockwise and Z(a, b) = arccos J^°' is the angle between the two vectors o, b G R 2 . 



j|a|| 

By o(e), e > 0, we will denote a function, which may depend also on other variables having the property that 

^ — > as e — > uniformly with respect to the other variables when they belong to a bounded set. 
Finally, let (,reM and let h(t, r) be the vector of R 2 given by 

h{t,r)=x Q {t)+r^§— (13) 
llM*) II 

Define the function (t, r) — > /(t, r) as follows 

J(t,r)=fi(T,0,ft(*,r)). (14) 

It is easily seen that, for any t G [0,T], the curve r — ► I(t,r) intersects the limit cycle xo at the point 
I(t,0) = x (t). 

The following theorem states a property similar to (9) in the case when the autonomous system (3) is perturbed 
by nonsmooth functions (f>. 

Theorem 1. Assume conditions (10). Assume that, for all sufficiently small e > 0, system (1) has a T -periodic 
solution x £ satisfying 

x s (t — 9o) — > x (t) as e — > 0, (15) 
for any t G [0, T], where 8 G [0, T]. Then there exist constants < Mi < M 2 , £0 > and r G (0, 1] such that 

eAfi|/i(0 o ,*)| < \\x e (t - e e (t)) -x {t)\\ +o{s) < sM 2 \f 1 (6 ,t)\ for any e G (0,e ) and any t G [0,T], (16) 

where 9 £ {t) — > 8 as e — > uniformly with respect to t G [0,T], and x e (t — # e (i)) G /(£, [— r ,r ]), i G [0, T]. 

To prove Theorem 1 we need some preliminary lemmas. 
Lemma 1. For any t G R we have 

im^Mt)z 1 {t))=\ {H0)M0)) I- ("0 

<yi(0),*i(0)> 



Proof. By Perron's lemma [32] (see also Demidovich ([6], Sec. Ill, §12) for any t G R we have 

(± (t) vi(«))*(*>(*) «i(*)) := 



(io(*),«o(t)> (io(*),«i(*)> I = / (io(0),«o(0)) (i (0),«i(0)> 
(yi(i)^oW) (yi(*),*i(t)> / \ (yi(0),z o (0)) (yi(0),«i(0)> 



Thus, in particular, (±o(0), zi(0)) = (io(7 1 ), zi(T)) . On the other hand Ao(0) = io(T) and zi(T) = p*zi(0), 
p* 7^ 1, thus (io(0), zi(0)) = 0. Analogously, since yi(T) — pyi(0), p ^ 1, and z (0) = Zo(T'), we have that 
(tfi(0),*b(0))=0. 

□ 
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Lemma 2. Under the assumptions of Theorem 1 there exist ro G (0, 1] and ao G [0,n/2) such that 



Z(I(t,r) - x (t), ±o(t) ) < a for any t G [0,T] and any r G [-r ,r ]. 



(18) 



Proof. Assume the contrary, hence there exist sequences {t n }neN C [0, T], t n — > to as n — ► oo, {r„} ne N C (0, 1], 
r„ — > as n — > oo such that 



Z(I(t n ,r n ) - Xo{t n ),x {t n )^) 7r/2 as n -> oo. 



We have 



I{t n ,r n ) - x (t n ) = fl (T,0,x {t n )+r n .. Z °^ n \ ± .. ) - x (t n ) = r n ^(T,0, x (t n )) ^4Htm + °( r ")" 

V F0(*n) 11/ ||M*n) || 



(19) 



(20) 



By Theorem 2.1 of [18] it follows that fy(T, 0, h(t, 0)) = Y(T, t) where Y(-, t) is the fundamental matrix for the 
system 

y(T)=^(x (T + t))y(r), (21) 



satisfying Y(0, t) = /, thus 0^(T, 0, 0))yi(t) = pyi(t). On the other hand from Lemma 1 we have 

Vi(t) II zo(t)^, 

therefore 0^(T, 0, 0))^o(i) ± = pzo(i) 1 " an d (20) can be rewritten as follows 

I(t„,r n ) - x Q (t„) = pr n ^r7TTT\\ +°( r «)- 



\\zo{tn)H 



(22) 



(23) 



Hence 

Z(I(t n ,r n ) - x (t„),x (t n ) ± ) = arccos 
Without loss of generality we may assume 



/ zoitn) 1 - ! o(r„ 
\ ll^o^n)^!! r n 




zoitn) 1 - o(r„) 


■llio^n)^!! 



and so 



(i o (0),z o (0)) = 1 



Z(I(t n ,r n ) - x a (t n ),x {t n )^) -> arccos - — 



(24) 



^(to^ll-llio^ll 

contradicting (19). Therefore there exist r G (0, 1] and a G [0,7r/2) satisfying (18). 



as n — > oo 



□ 



We can now prove the following. 

Lemma 3. Under the assumptions of Theorem 1 there exists so > such that for any e G (0, eo) and any 
t G [0, T] we have 

x £ (t-6 e (t)) G I(t, [-ro.ro]) 

w/iere e (i) = 6 - A e (i), A e (t) G [i - f ,* + f] and A e (i) — > as e — ► uniformly in t G [0,T]. Moreover, 
there exists M > swc/i i/iai 

||x £ (t - e (t)) - x (i)|| < eM /or any £ G (0, e ) and any £ G [0, T]. (25) 
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Proof. First of all observe that ro > given by Lemma 2 can be chosen to satisfy 



/(t,[-r ,r ])na;o(M) = {a;o(t)}, for any t G [0, T] . (26) 

From (18) of Lemma 2 and (15) we have that there exists £o > such that I(t, [— ro,ro]) H x £ ([0,T]) / for 

any e G (0,£o) and an y * G [0, T\. Hence, for any e G (0,£o) an d * £ [0, T] there exists A e (i) G [t — + -j] 
such that 

^(t-^o + M*)) [-ro.ro]). (27) 

We claim that 

A £ (t) -> as e -> (28) 

uniformly with respect to t G [0, T]. In fact, assume the contrary, thus there exist sequences {e n } ne tq C (0,e ), 
e n — > as n — > oo, and {i„}„ e N, tn — > *o G [0,T] as n — > oo, such that A £(l (t„) — > A ^ and x £n (t n — 8 a + 
A £n (i„)) G J(t n , [-r ,r ]). Since from (26) we have 7(t n , [-r ,r ]) nx ([0,T]) = {x (i„)} then 

x £n (t n -6q + A £n (t n )) -> x (t ) as n -> oo. (29) 

Applying (15) we have 

^„ (^ - #o + A £ „ (t„)) -» x (A + t ). (30) 
From (29) and (30) we conclude that 

zo(to) -x (Ao + i ), (31) 

where A G [to — 5, to + -j] , since T is the smallest period of xo it follows from (31) that A = 0, which is a 
contradiction. 

Pick any r G [0, T], in what follows we show that the shifts t — » A £ (i) have the property that the convergence 
of x £ (r + i — # £ (t)) to x (r + t) is of order £ > 0, where E (t) = 9 — A £ (t), and thus the claim of Lemma 3 is 
proved. 

For this consider the change of variables u £ (r,t) = Cl(0, r, x £ (r + t — 6 £ (t))) in system (1). It is clear that 
x £ (t + t — £ {t)) = fi(r, 0, v £ (t, tj) and so 

x e (r + t- e e (t)) = V(fi(r, 0, v e {j, t)) + Q' € ( T , 0, i/ e (r, t))(i/ e );(r, t). (32) 

On the other hand from (1) we have 

x £ (t + t - e e (t)) = i/>(Sl( T , 0, i/ e (r, i))) + £0(t + i - e (i), fi( r , 0, i/ e (r, t)), e). (33) 

From (32) and (33) it follows 

(i/ e );(r, t) = £ (fi^r, 0, i/ e (r, t))) _1 0(t + t - 9 £ (t), fi( T) 0, i/ e (r, t)), e) 

and since 

u £ (0,t) =x £ (t- 6 £ {t)) =x £ (T + t-e £ (t)) = Sl(T,0,v e (T,t)) (34) 

we finally obtain 

u e (T,t)=n(T,0,v e (T,t))+e f (n£(a,0,i/ e (s,t))) _1 <^(a + t - 6 e (t),Sl(a,0,u e (s,t)),e)d8. (35) 
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Since v e (t, t) — > 0(0, r, xo(t + t)) = xo(i) as e — > we can write ^ e (r, t) in the following form 

u E ( T ,t) =x (t) +£H e {T,t). (36) 
Subtract Xo(t) from both sides of (35) obtaining 

en e (T,t) = eQ' i (T,0,xo(t))p E (T,t) + o(ep E (T,t)) + 

+ e [ (nl(s,0,i/ e (a,t))) _1 <j){s + t - 9 E (t),Sl(s,0,"e(s,t)),e)ds. (37) 
Jo 

Since x e (t — e (t)) £ I(t, [— r , r ]) then from (14) there exists r e (i) £ [— r ,r ] such that x e (t — e {t)) = 
Q.(T,0,h(t,r E (t))) and by (13) we get 

eii e (T,t) = v e (T,t)-xo(t) = Sl(p,T,x e (t-6 e (t)))-xo(t) = 

zoit) 1 - 



= fi(0, T, fi(T, 0, h(t, r e (t)))) - x (t) = h(t, r e (t)) - x (t) - r e (t). 



Therefore fi e (T,t) || ^(i) -1 and by (22) we can rewrite (37) as follows 

eHe(T,t) = epp E (T,t)+o(ep E (T,t))+s f 0, v e (s,t)))~ 1 <p(s + 1 - e e (t),Sl(s, 0, v e (s, t)),e)ds. (38) 

Jo 

We now prove that the functions (r, t) — ► /i e (r, i) are uniformly bounded with respect to e £ (0, e ). For this we 
argue by contradiction, therefore there exist sequences {£ n } n eN C (0,1), e n — > as n — > oo, {t„}„ £ n C [0,T], 
t„ — > To as n — > oo and {t„}„ e N C [0, T], — > to as n — > oo, such that ||/ie n (r n ,t n )|| — > oo as n — > oo, so 
1 1 Me (• ) tn)||c — > 00 as n — > oo, where || • ||c is the usual sup-norm of C([0, T], M 2 ). Let g n (r) = f^lnjll ^ ; 
then from (38) we have 

/ \ /rri\ I o(e n p En (T,t n )) 

qn(r) = pq n {T) + + 



1 



(n^(s, 0, (s, t n ))) 1 0(s + t„ - £ „ (t n ), Sl(a, 0, ^ £ „ (s, t n )), e n )ds. (39) 



||Me n (-)*n)||c 7o 

By definition the set of continuous functions A = {q n , n £ N}, is bounded and, as it is easy to see from (39), 
A is also equicontinuous. Therefore, by the Ascoli-Arzela Theorem, see e.g. ([4], Theorem 2.3), we may assume 
without loss of generality that the sequence {g n }neN is converging. Let qo = limn-nx, q n , from (39) we may 
conclude that 

qo(r) = pqo(T). (40) 

By (40) it follows that q is a constant function, thus being p ^ 1 we have qo = 0. On the other hand, by the 
definition of q n , we have that || (folic = 1- This contradiction shows the uniform boundedness of the functions 
p E with respect to e £ (0, £o). On the other hand from (34) and (36)wc have that 

x e (t-9 e (t))-x (t)=ene(0,t), (41) 

and thus the proof is complete. 

□ 

Proof of Theorem 1. We have to prove (16) with t — > e {t) as given in Lemma 3. For this, by Lemma 1, we 
can represent x e (r + t — 6 £ {t)) — x (t + t) as follows 

x e (t + t - 6 E (t)) -x (r + t) = ea E {T,t)x (T + t) + eb E {T,t)yi(T + t), (42) 
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where 

ea £ (r,t) = (z (t + t), x £ (t + t - £ (t)) -x (T + t)) and 

sb £ (T, t) = ( Zi (t + t),x e (r + t- 6 £ {t)) -x {t + t)) . (43) 

By Lemma 1 we have that (xa(t), z\{i)) = 0, for any t e [0, T], and so xo(t) x = k - ' )) z±(t), where k = +1 

\\ z i\t)\\ 

or k = — 1. Therefore 

(±o(r + t)^,x £ (r + t - 9 £ (t)) - x (r + t)) = eb £ (r, t)k ^ o(t + t) J~f . (44) 

\\Zl(T+t)\\ 

We aim now at providing an explicit form for (44) by looking for a suitable formula for the function (r, t) — > 
b £ (r,t). To do this we substract (3) where x(r) is replaced by x (r + t) from (1) where x(t) is replaced by 
x £ (t + t — 9 £ (t)) to obtain 

x £ (t + t - 6 £ (t)) - x (r + 1) = iI/(x (t + t))(x £ (r + t - 9 £ (t)) - x (r + t)) + 

+ s4>(t + t - 6 £ {t), x £ {t + t - 6 £ {t)), e) + o(x £ (t + 1 - 6 £ {t)) - x (t + t)). (45) 

By substituting (42) into (45) and taking into account that 

ea £ {r, t)ip'(xo(T + t))xo(r + t) = ea £ (r, t)xo(r + t) and 
eb £ {r, t)ij}'{x {T + t)) yi ( T +t) = eb £ (r, t)yi{r + t) 

ex (T + t)(a £ )' T (T, t) + e yi (T + t){b e )' T (r, t) = £0(r + * - 9 £ (t),x £ (T + t - 6 £ {t))) + 

+ o(x £ (t + t - 9 £ (t)) - x (t + *)), 

e(b £ y T (T,t) = e(z 1 {r + t),cj ) (T + t-9 £ (t),x £ (T + t-e £ (t)))) + 

+ ( Zl {T + t),o(x £ {T + t-e £ (t))-X {T + t))}. (46) 

Moreover, since Zi(t) = p*zi(r — T), from (43) it follows that 

b £ (T,t)=p*b £ (T-T,t). (47) 
System (46)-(47) has a unique solution which, as it is easy to verify, is given by the formula 
b £ (T,t) = J^jJ t (zi{s + t),(f)(s + t-e £ {t),x £ {s + t-e e (t)),s))ds + 

P* - 1 Jr-T \ £ / 

By substituting this formula into (44) we obtain 



we have 



and so 



= ek 



{x {T + t)^,X £ {T + t-e £ {t))-X {T + t)) = (48) 



|*i(T + t)||(p.-l) 



f (z!(s + t), 4>{s + t - 9 £ (t),x £ (s + t - 9 £ (t)),e)) ds + 

Jt-T 



IMr + ^i \p. r ( Zl{s + t) ^ + t-e £m - M s + t)y 

\\zi(T + t)\\(p*-l) J T _ T \ e 



where °^ x ^ s + t — 6 J£>} x o( s + *)) ^ q as e ^ uniformly in s G [-T, T] in virtue of (25). On the other 



hand 



{x {T + t)^,x e {T + t-6 £ (t)) -x (r + t)) = 
= \\x {T + t)^\\ \\x £ {T + t-6 e {t))-x (T + t)\\ cosZ(i (T + t)- L ,a; e (T + t-fl e (t)) -x (T + t)) 



and by taking into account (18) of Lemma 2, (48) and the fact that ||i;o(t)" L || ^ for any t G [0, T] we get 

rO 

-T 



\x £ {t-6 £ {t)) -x (t)\\ = eg{t) J ^{z 1 {s + t),4>{s + t-6 s {t),Xe{s + t-6 e {t)),s))ds + 



f° I o(x £ (s + t-9 £ (t))-x (s + t)) \ , 
£9{t) (zi{s + t), )ds, 



where 



9(t) 



\\zi(t)\\(P* - 1) cosZ(i o (t) x ,a; e (t-0 e (t))-a;o(t)) 
is a continuous function on [0, T] with g(i) ^ for any t G [0, T], Therefore, 

||a; e (t-0 e (t))-a;o(t)|| = eg{t) j { Zl {s + t),<p(s + 1 - 6 £ {t),x £ {s + t - 6 e (t)),e)) ds + o(e). 



(49) 



On the other hand from Lemma 3 we have that A e (i) — > uniformly in t G [0, T], thus we can rewrite (49) as 
follows 

\\x s {t - 6 £ {t)) - x {t)\\ = eg{t) J ( Zl (s + t), <f>{s + t- 6 , x {s + t), 0)) ds + o(e), 
introducing the change of variable s + t = u in the integral we finally get 

\\x £ {t-6 £ {t)) -x (t)\\ =eg(t) f (z^u), <f>{u - 6 , x e (u - 6 ), 0)) du + o{e) 

Jt-T 

from which (16) can be directly derived recalling that, by Lemma 3, x e (t — £ (tj) G I(t, [— ro, ro]) for any 
e G (0, s ) and any t G [0, T). 

□ 

As a straightforward consequence of Theorem 1 we have the following result. 

Corollary 1. Assume all the conditions of Theorem 1, then for every t G [0,T] such that 

fi(0 ,t)=0 

we have 

\\x £ (t - e £ (t)) -x (t)\\ =o(e), 
where £ {t) — > #o as e ^ and x e (i — 6* e (i)) G J(i, [— ro,ro]). 
Next result is also a consequence of Theorem 1. 

Corollary 2. Assume all the conditions of Theorem 1. Moreover, assume that 

A (0 O ,*) 7^0 for any t G [0,T]. (50) 
TTien i/iere ezzsis ei > swc/i that x £ (s) ^ xo(t) for any s,t G [0, T] and any e G (0,ei). 
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Proof. Let £o > given by Theorem 1. From (50) we can choose E\ G (0,£o) in such a way that, for any 
e G (0, £i), we have both 

o(e) < eMi|/i(0 o ,t)|, for any * G [0,T], (51) 

and the validity of (16). Moreover, E\ can be also chosen in such a way that there exists So > such that 
the curve t — > a; e (r) intersects /(t, [— #o,£o]) a ^ om y 0110 P om t for any e G (0,£i) and t G [0,T]. Such a choice 
is possible, in fact, since x s {t — 9o) — > io( T ) as e uniformly with respect to r G [0,T] and the curve 
r — > J(t, r) intersects the limit cycle xo transversally at r = for any t G [0, T], then there exists 5o > such 
that x e ([t ~9o — So,t — 9o + So]) and I(t, [—So, So]) have only one common point for any t G [0, T] and sufficiently 
small e > 0. On the other hand r — > £ £ (t) cannot intersect /(t, [— <5 ,<5o]) for r G [i — #o — ^, t — #o — So] U 
[t — #o + <5o> t — #o + 5] and e > sufficiently small, otherwise there would exist sequences {e„}„eN, e n — > as 
n -> oo, {t„}„ eN , t n —> t Q G [0,T] asn^oo, {r n }„ eN , r„ G [t„-#o - *n - #o - <>o] U [t„ - 9 + S , t n - 9 + f ], 
t„ -> t G [t - #o - f , *o - #o - $o] U [t - #o + S , t - 9 + f ] as n -> oo such that x £n (r„) G J(t„, [-<5 , <5 ]), 
thus x (t + Oo) = Xo(to), with r + #o 7^ to an d |t + 6* — to| < T, which contradicts the fact that T > is the 
smallest period of x - 

To conclude the proof assume now, by contradiction, that there exist e G (0,£i) and s, i G [0,T] such that 
£e(s) = so(t)- Since t — > x £ {t) intersects I(t, [— So, So]) at only one point then Theorem 1 implies that s = 
t — 9i(t). In conclusion, from (16) we have eM\\f\{9o, t)\ < o(e) contradicting (51). 

□ 

The following result is crucial for the proof of our existence result Theorem 3, but it can be also considered as 
an independent contribution to the coincidence degree theory. 

Theorem 2. Assume conditions (10). Fore > 0, let G £ : C([0,T],M 2 ) -> C([0,T],R 2 ) be the operator defined 
by t t 

(G e x)(t)=x(T)+ [ i>(x{T))dT + e [ cP(r,x(T),e))dT, te[0,T]. 
Jo Jo 

Let W Ua = {x G C([0,T],R 2 ) : x(t) G U , for any t G [0,T]} . Assume that 

(io(0),z o (0)} = (y 1 (0),z 1 (0))-l. (52) 

Finally, assume that for every 9 G [0, T] such that fo(&o) = we have 

fi(0o,s + 9 )^0, for any s G [0,T]. (53) 

Then, for all e > sufficiently small, I — G £ : C([0, T],R 2 ) — ► C([0,T],R 2 ) is not degenerate on the boundary 
of Wjj ■ Furthermore, there exists a continuous vector field F : R 2 — > R 2 such that 

d(I - G e ,W Uo ) = d B (F,U ), 

where F(x (9)) = fo{9)xo{9) + 9) Vl (9) for any 9 G [0, T]. 
Some remarks are in order. 

Remark 1. As already observed condition (52) does not affect the generality of Theorem 2. 

Remark 2. In Theorem 2 we could replace dB{F,Uo) by k ■ ind(xo, F), where k = +1 or k = —1 according 
with the orientation of the limit cycle xq- Precisely, k = +1 if the set Uq is on the left side when one follows 
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dUo according to the parameterization xo(t) with t increasing from to T, and k = — 1 in the opposite case. 
Moreover, u\d(x ,F) is the Poincare index of the trajectory xo with respect to the vector field F, namely the 
total variation of an angle function of the vector F(xo(t)) when t increases from to T, see Lefschetz ([21], 

Ch. ix, § 4). 

Remark 3. The Jordan theorem, see Lefschetz ([21], Theorem 4-V> ensures that the interior Uq of xq does 
exist and it is an open set. 

To prove Theorem 2 we need the following preliminary lemma. 
Lemma 4. For any s e [0, T], let 

F„(Z)= f fi£(0, t,O(t,0^))^(t,O(t, 0,0, 0)dr, for any £ e K 2 . (54) 

Js-T 



Then 



(zo{6),F s {x (6))) = fo(0) for any s,0G[O,T], 

(z 1 {0),F s {x o (0))) = f 1 {0,s + 0) for any s,0G[O,T}. (55) 



In particular, 



F,(x (9))= 1 — f o {9)xo{0) + - 1 f 1 (0, S + 6)y 1 (6) for any s,0,t e [0,T]. (56) 

{x (t),z {t)} {yi(t),zi(t)) 

Proof. It can be shown, see Krasnosel'skii ([18], Theorem 2.1), that fl'^(t, 0, xq(6)) = Y(t,0), where Y(t,0) is 
the fundamental matrix of the system 

y(t)=*['(x (t + 9))y(t) (57) 
satisfying Y(O,0) = I and since fi£(0, t, fl(t, 0, x (9))) ■ n' s (t,O,x o (0)) = I we have 



F s (x (9)) = f Y- 1 (T,9)<j>(T,x (T + 9),0)dT. 

Js-T 



(58) 



Let us now show that 

F^M) = Y(9,0)Y- 1 (t + 9,0). (59) 

In fact, it is easy to see that Y(t + 9,0) is a fundamental matrix for system (57) and so Y(t + 9,0)Y~ 1 (9,0) 
is also a fundamental matrix for (57), moreover we have that Y(t + 9,0)Y~ 1 (9,0) = I at t = 0. Therefore 
Y(t + 9,0)Y- 1 (9,0) = Y(t,9) which is equivalent to (59). 

By substituting (59) into (58) and by the change of variable r + 9 = t in the integral of (58) we obtain 
F s {x o (0))=Y(0,O) [ Y~\t + 9,0)^(t,x (t + 9),0)dr = Y(0,O) / + Y' 1 (t,0)^(t - 9, x o (t),0)dt. 

Js-T Js-T+6 

Let Z(t) be the fundamental matrix of system (6) given by Z(t) = Zo(t)Z ( ^" 1 (0), where Z^(t) — (zo(t) z\{t)), 
t G [0,T]. Since y _1 (i,0) = Z*(t), see Perron [32] and Demidovich ([6], Sec. Ill, §12), then we have 

F s (x a (9))=Y(9,0) [ S+e Y-\t,^{t ~ 9,x {r),0)dr = {Zmy 1 Z*(t)<[(t - 9, x o (r),0)dT. 

Js-T+e Js-T+e 

Let 

rs+e 

A(s,0)= Z*(r)cP(T-9,xo(T),0)dT, 
Js+e-T 



is+e-T 

12 



we have 

' 'zi(0)] 2 / I [zi(0)]i [*i(*)] 2 



2 



A( S , 



k(0)] 2 j ' dctZo(0) ^ [zo(e)]i 
1 // \ / [zi(e)] 2 [A(M)]i-NW]2[A(M)] 2 

{[ Zl {e)Wi{o)HHsMi - [zi{6)Uz {eM&{s,6)]2 - 



detZo(fl) 

- ^W] 2 [^ 1 (^)] 1 [A( S ^)] 1 + ^W] 2 ^oW] 1 [A( s ,e)] 2 } = 
= {[A.(fl)]i[«i(9)] 2 - [^oWM^^J [A( S , 6)} i+1 = 

1 /- s+e 

= A +7 M detZ °W / (^W.^- Mo(t),0)) dr. 
detZ (6>) A+o-t 

For i = 0, 1, we obtain (55). Furthermore, from Lemma 1 and (52) we have that 

F s (x (6)) = 1 — f Q (O)x Q {0) + 1 + %i(0) for any G [0,T] and any i G [0,T]. 

(io(i),2b(<)) (2/iW)^iW) 

□ 

Proof of Theorem 2. 

Let 77(4, s, £) be the solution of the system 

q(t) = i/(ii(t, 0, 0)q(t) + n(t, 0, 0) (60) 

satisfying ?7(s, s,£) = whenever £ G R 2 . It can be shown, see ([16], Lemma 2), that 

F s (0=v(T,s,®-r 1 (0,s,Q. (61) 

Therefore, from (52), (53) and (56) we have that 

r?(T, s, £) - 77(0, s, 0^0 for any £ G <9f7 and any s G [0, T], (62) 

and by applying ([15], Theorem 2) we obtain the existence of an £0 > such that 

d(I-G e ,W Uo ) = d B (v{T,0,-),Uo) for any eG(0,£ ), 

where W Uo = {x G C([0,T],R 2 ) : fi(0,t,a:(t)) G J7 for any t G [0,T]}. Since as it is easy to sec W Ua = W Ua , by 
taking into account (56) and (61) we end the proof by defining F(£) = i*b(£), £ G R 2 . 

□ 

The following existence theorem is the main result of the paper. 

Theorem 3. Assume (10). Assume that for every zero 9 G [0,T] of the bifurcation function f we have 

/1 (00, *) 7*0 for any t G [0, T]. (63) 

Let F G C(R 2 ,R 2 ) be a vector field such that on the boundary of Uo it has the form F(x (9j) — fo(6)xo{6) + 
f 1 (9,9)y 1 (9) for any 8 G [0,T]. Assume 

d B (F,U )^l. (64) 
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Then there exists £o > such that for every e S (0,£o) system (1) has at least two T -periodic solutions x\. e and 
X2, E satisfying 

%i,e(t — 9i) — > x n (t) as e — > 0, i = 1,2, (65) 
w/iere #i, # 2 G [0, T]. Moreover, we have that xi, £ (t) E Uq and x 2tE (t) & Uo, for any t E [0, T] and any e E (0, £o). 

Proof. Denote by WsidUo) the (5-neighborhood of the boundary 8Uq of the set Uo- Let Ui t s = Uo\Ws(dUo) 
and U 2 ^ = Uo U Ws(dUo), thus the set C/i ;< 5 tends to Uo from inside as 6 — > 0, while [^2.5 tends to C/o from 
outside as <S — > 0. Since the limit cycle xo is isolated then there exists S > such that 

Gq(x) 7^ x for any x E dWu l s U dWu 2 S and any Se(0,S ]. (66) 

Moreover, being T > 0, we can choose <5 > in such a way that 

V>(£) + for any £ G U 9J7 2 ,5 and any (5 e [0, <5 ]. (67) 

From (67) we get 

dB(i>, Ui,s ) = d B {ip, U 2 ,s ) = d B (ip, U ). 

Since Uo is the interior of the limit cycle Xo of system (3) by Poincare theorem, see Lefschetz ([21], Theorem 11.1) 
or Krasnosel'skii et al. ([19], Theorem 2.3) we have d B {ip, Uo) = 1 and so 

ds(V>, Uija) = d B (ip, U 2 ,s ) = 1. 

In virtue of (66) and the fact that Wu n K 2 = U, ([5], Corollary 1) applies to conclude that 

d(I - G , W Ui Sq ) = d B (ip,U liSo ) and d(I - G , Wjj 2Sq ) = d B (ip, U 2 ,s ), 

hence 

d(I - G ,W Ult J = 1 and d(I - G , W v ^ o ) = 1. 
Therefore, there exists £0 > such that 

d(I - G £ ,W UliSo ) = d(I - G £ ,W U2i J = 1 for any £G (0,£ ). (68) 

Since by the definition of z\ we have that z\(t + T) = p*zi(t) for any t E [0, T], then for any t E [0, T] it is easily 
seen that f\{9, t + T) = p*fi(9, t), whenever 9 E [0, T], and thus from (63) we have also that /i(#o, t + 9 ) ^ 
for any t E [0, T}. Therefore, all the conditions of Theorem 2 are satisfied and we can take £0 > sufficiently 
small to have 

d(I-G e ,W Uo ) =d B (F,U ) for any ££ (0,£ ). (69) 
By (64), (68) and (69) we conclude that for any e E (0, £0) there exist 

xi, e €W Uo \W UliSo , and x 2 , e E Wjj 2Sq \W Uo (70) 

such that G E (xi^ e ) = x\ tS and G £ (x 2 , £ ) = x 2 . e - From (70) we have that for any e E (0,£o) there exist points 
h, s , t2, £ G [0, T] such that x ly£ (t 1:£ ) E U \U 1A> and x 2y£ (t 2 , £ ) E U 2A> \U . Thus xi. s {t) —> dU and x 2 , £ (t) —>■ dU , 
for any t E [0,T], as e — > 0, otherwise there would exist a T-periodic solution x* to system (3) and a point 
t» E [0,T] such that either x»(t») £ t^oWiA or x *{t*) G U 2 ,s \Uo contradicting (66). Therefore, see ([25], 
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Theorem p. 287) or ([23], Lemma 2), for every i £ {1,2} there exists 0% £ [0, T] satisfying (65). The fact that 
Xi t e(t) £ Uo and x 2 , e {t) Uq for any t £ [0, T] and e > sufficiently small follows from Corollary 2 and so the 
proof is complete. 

□ 

Remark 4. From the proof of Theorem 3 it results that d(I — G £ ,Wu \Wu lt g ) and d(I — G e , Wu 2 ^ \Wu ) 
are different from zero for e £ (0, £o). This fact can be used to obtain stability properties of solutions Xi iE and 
x 2 , e in the case when further information on the number of T -periodic solutions to (1) belonging to the sets 
WuoXWux^a an d Wu 2 ,6 \Wu ) are available, see Ortega [31]. 

3. An example. 

In this section we always assume that condition (A ) is satisfied. The well known formula by Poincare, see 
Krasnosclskii ct. al. ([19], formula 1.16) states that 

md^F)--^ lamhTWm ^ 

where a(9) — F(xo{0)), 9 £ [0, T]. The relationship between ind(xo, F) and d B (F, Uq) was discussed in Remark 2. 
In this section we show how the representation F(xo(8)) = fo(6)xo(6) + fi(6, 9)yi(0),0 £ [0, T], of the function F 
on dUo permits a simpler calculation of ds{F, Uq). For this, we consider the case when <f>(t,£) = —<j)(t + T/2,£), 
which includes, in particular, the classes of perturbations <f>(t, x) = smt ■ (f>i(x) and <j){t, x) = cost • (f>i(x), where 
<t>i £ C(R 2 ,R 2 ). We can prove the following result. 

Proposition 1. Let F £ C(R 2 ,M 2 ) be a vector field such that F(x o (0)) = f o (0)x o {0) + fi{0,T)x^(0), 8 £ [0,T]. 
Assume that 

(xo(0),^(0)> = (i/i(0),«i(0)> = 1, (71) 

fo(O) = -fo(0 + T/2) for any 9 £ [0,T], (72) 

f 1 (0,T) = -f 1 (0 + T/2,T) for any £ [0,T]. (73) 

Moreover, assume that there exists an unique 9q £ [0, T/2) such that fo(&o) = 0. Finally, assume that the 
function f is strictly monotone at the point 9q and that 

/i(0o,T)^O. (74) 

Then either d B (F, U ) = or d B (F, U ) = 2. 

The proof of the proposition is based on the following technical lemma. 

Lemma 5. Let U C M 2 be an open set whose boundary dU is a Jordan curve q : [0, T] — > M 2 , with q(0) = q(T). 
Let F : M 2 — > K 2 be a continuous vector field such that F(£) ^ for every £ £ dU. Assume that for a continuous 
function z : [0, T] — > R 2 , z(0) = z(T), the following conditions hold: 

1) (z(9),q(9)) ^ for every 9 £ [0,T], 

2) the function f(9) = (z(0),F(q{0))^ has exactly two zeros 9i,6 2 £ [0,T), 

3) the function f is strictly monotone at 6\ and 2 , 

4) sign («(fli) x , F(g(&i))) - -sign (z(6 2 ) L , F(q(9 2 ))) . 
Then either d B (F, U) = or d B {F, U) = 2. 
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Proof. Assume that the parametrization q is positive, namely the set U is on the left side if one follows 
dU according to the orientation given by q when t increases from to T, otherwise we consider the opposite 
parametrization q(9) = q{—6). For any t £ [0,T] we denote by Q(t) the angle (in radians) between the vectors 
q(0) and q(t) calculated in the counter-clockwise direction. Clearly @(t) is a multi-valued function of t. Let 
Tg(t) be the single- valued branch of 9(t) such that T q (0) = and let Q : dU — > M 2 be the vector field defined 
by Q(q(t)) := q(t), whenever t G [0,T], hence T q (t) = T Qoq (t). Following ([19], §1.2) the function t -> T q (t) 
is called the angle function of the vector field Q on the curve q . Analogously, considering the angle between 
F(q(0)) and F(q(t)), we can define the angle function Fp oq (t) of the vector field F on the curve q. By the 
definition of the rotation number for planar vector fields on the boundary of simply-connected sets, see ([19], 
§ 1.3, formula 1.11) we have 

dB(F,U) = ^[Tp oq (T)-Tp oq (0)}. (75) 

Therefore, in order to prove the lemma we must calculate the right hand side of (75). For this, denote by hi, hi £ 
[0, 27r) the angle between the vectors hi and hi in the counter-clockwise direction, that is hi, hi + hi, hi = 2ir. 
Observe that 

r i?o 9 W - r «W = r q ^ oq (0) mTF{q{0)), (76) 
where ^ q p oq (0) is the single valued branch of the multi-valued angle between q(9) and F(q(9)) such that 
r q ^ oq (0) = q(0),F(q(0)). 

To calculate T q Foq (9) we introduce the function Z : M 2 x M 2 — > [—n, n] as follows 



hi, hi as hi, h 2 G [0,7r], 

hi,hi — 2iT as /ii, /i 2 G (7T, 2-jt] 



By condition 3) we have that ind(#i,/) = +1 or ind(0j,/) = —1 according to whether / is increasing or 
decreasing at 6i, i — 1,2. 

Up to a shift in time, since 9 2 — &i < T, we may assume that the zeros 9i,6i of f{6) = (z{6), F(q(9))^ 
belong to the interval (0, T). 

Assume that (z(6),q(6)) > for every 9 £ [0, T], otherwise we consider z(9) = z{—9) instead of z(9). A 
possible way to write explicity the function T q Foq {9) is the following 



r i,i?o,(0) = { 



Z(z(9),q(9)) + ^(sign (z(9), F(q(6))) z(9),F{q{9))) as 9 £ [0, 0i), 
^{z{6),q{e)) + /(sign (z{6), F(q(9))) z{9),F{q{9)))+ 

+7rind(# 1 ,/)sign(z(0 1 )\F(< Z (# 1 ))) as 9 £ {61,61), 
^{z{6),q{6)) + Z(sign (z{6), F(q(9))) z{9),F{q{9)))+ 
+7rind(fli, /)sign ^(flO-L, + 

+7rind(0 2 ,/)sign(z(0 2 )\.F(g(0 2 ))\ as 9 £ {6 2 ,T\. 



It is easy to see that the above representation of the function 9 — > T. p oq {9) can be extend to 9i and 9 2 by 
continuity. Since 

- /(sign /z(0), 
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are T-periodic functions from (75)-(76), taking into account that 

d B (Q,U) = ^-[T i (T)-T 4 (0)] = l, (77) 

(sec e.g. [19], Theorem 2.4), we have 

d B {F,U) = l+ l - [ind(0 1 ,/)sign(z(0 1 ) ± ,F(^ 1 )))+ind(0 2 ,/)sign(z(0 2 ) ± ,^( (Z (0 2 )))] . (78) 
Since the function / is T-periodic then 

ind(0i,/) = -ind(0 2 ,/) (79) 
By assumption 4) and (79) the claim can be easily derived from (78). 

□ 

Proof of Proposition 1. 



xp(t) 
\\x (t)\ 

in Lemma 5. Let us now show that all the conditions of Lemma 5 hold. In fact, we have that (io(t), z(t)) = 1 



Let U = Uq, q(t) = xo(t), z(t) = — — , t G [0, T], thus the function /o turns out to be the function / defined 

/ 1 1 z 



for any t G [0,T] and so condition 1) is satisfied. Our assumptions imply that the function fo has only two zeros 
9\ = 9q and 02 = 0q + T/2 in the interval [0, T] and it is strictly monotone at these points, thus conditions 2) 
and 3) of Lemma 5 are also satisfied. Finally, (^(6)^-, F(xq(0))^ = fi(6,T) and so (73) implies condition 4) of 
Lemma 5. Hence the proof is complete. 

□ 

By combining Theorem 3 and Proposition 1 we obtain the following result. 
Corollary 3. Assume conditions (10) and assume that 

<f)(t, £) = -<j>(t + T/2, £) for any t G [0, T] and any (el 2 . 

Moreover, assume that there exists a unique 9q G [0, T/2) such that /o(#o) = 0. Finally, assume that the function 
fo is strictly monotone at the point 9q and 

fi(0o,t)?O for any t G [0,T]. (80) 

Then there exists Sq > such that for every e G (0,£o) system (1) has at least two T-periodic solutions and 
X2. 6 satisfying 

Xi, e (t - 0i) — > x (t) as e — *■ 0, % = 1, 2, 

where 81,62 G {0q,9q + T/2} . Furthermore, we have that X\, s (t) G Uq and X2, s {t) £ Uo, for every t G [0, T] and 
e G (0,e ). 

Proof of Corollary 3. 

To apply Theorem 3 we only have to verify condition (64) . For this we will make use of Proposition 1 . Without 
loss of generality we can assume that 

(l/i(0),x- L (0))>0. (81) 

We claim that, under the conditions of Corollary 3, the vector field F of Proposition 1 is homotopic on dUo to 
the vector field F of Theorem 3. To prove the claim we show that the following homotopy joining F and F 

D x {x {e)) = f {6)x a {9) + h{6,\T + {I - X)6){\x a {e)^ + {I - X) yi {9)), with 9 G [0,T] and A G [0,1], 
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is admissible. Assume the contrary, therefore there exist Ao G [0, 1] and 9o G [0, T] such that 
fo(Oo)x o (0o) + /i (0o, A T+ (1 - Xo)e o )(X o x o (0 o ) ± + (1 - Ao)l/i(0o)) = °- 
By condition (81) we have that the vectors xq(0o) an d \oXo{0o)~ L + (1 — Ao)j/i(0o) are linearly independent thus 

/o(«o) = and f 1 (e ,\ T+(l-\ )6 ) = 
contradicting assumption (80). Hence we have proved that 

d B (F,U ) = d B (F,U ). 

Applying Proposition 1 we obtain that 

d B (F,U ) G {0,2}, 

namely assumption (64) of Theorem 3 is satisfied and the conclusion of the corollary follows from Theorem 3. 

□ 

At the end of the paper we would like to stress that all the functions j/i, z and z\ can be easily determined 
both analytically and numerically once the limit cycle x is known. We give in the following a sketch of both 
approaches. 

1) The analytical approach. 

Since xq is one of the two eigenfunctions of system (4) then by using well known formulas, see e.g. Pontrjagin 
([33], p. 138), the dimension of the system (4) can be decreased by 1, thus the obtained one-dimensional system 
can be easily solved to determine y\ . Furthermore, by Lemma 1 the eigenfunctions zq and z\ can be determined 
by the formula 

( Zo (t) Zl {t)) = {{x Q {t) yi (t)rr 1 . 

2) A direct numerical approach. 

From Lemma 1 we have (xo(0), zi(0)) = 0, therefore as initial condition we may take zi(0) = xo(O) 1 - and then 
z\ can be obtained by a numerical computation. By the definition of z\ there exists a T-periodic function 
a G C(M,K 2 ) such that z\{t) = a(t)e p **. Assume, that < 0. Let us fix an arbitrary vector £ G K 2 , which 
is linearly dependent with zi(0) and calculate the solution z of system (6) satisfying z(0) = £ on the interval 
[0, kT] where k £ N. It turns out that larger is k better accuracy is obtained. Observe that z can be represent 

by 

z(t)=aa(t)e p * t + z (t), (82) 
where z is an eigenfunction of (6), and since e p *' — > as t — > +oo then for given k e N we may take 

z (t) = z(t+(k - 1)T) for any t G [0,T]. 

For the case p* > one should make the change of variables z(t) = z(—t) for any t G M. in (6), to calculate zq on 
[—kT, 0] (for this it is necessary to expand, the function z\ on the interval [—kT, 0]) and then put z (t) = z\{—t) 
for any t G [0,T]. Once zo is calculated with the desirable accuracy the function y\ can be determined as the 
solution of (4) with initial condition yi(0) = Zo(Q) ± . 
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